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Abstract 

Based on a chiral approach to nuclear matter, we calculate the quartic term in the expansion 
of the equation of state of isospin-asymmetric nuclear matter. The contributions to the quartic 
isospin asymmetry energy A 4 ^{kf) arising from Ivr-exchange and chiral 27r-exchange in nuclear 
matter are calculated analytically together with three-body terms involving virtual A(1232)- 
isobars. From these interaction terms one obtains at saturation density po = 0.16 fm“^ the 
value A 4 (A:/o) = 1.5 MeV, more than three times as large as the kinetic energy part. Moreover, 
iterated Ivr-exchange exhibits components for which the fourth derivative with the respect to 
the isospin asymmetry parameter 6 becomes singular at <5 = 0. The genuine presence of a 
non-analytical term (5^1n |5| in the expansion of the energy per particle of isospin-asymmetric 
nuclear matter is demonstrated by evaluating a s-wave contact interaction at second order. 


PACS: 21.60.Jz, 21.65.+f, 24.10.Cn. 


1 Introduction and summary 

The determination of the equation of state of isospin-asymmetric nuclear matter has been a long¬ 
standing goal shared by both nuclear physics and astrophysics [1] . Usually one assumes a parabolic 
form for the energy per nucleon at zero temperature, E^s^pp, pn) = E{p) + A. 2 {p) where 

P = Pp + Pn is the total nucleon density and 5 = — Pp)/p is the isospin asymmetry related to 

unequal proton and neutron densities Pp ^ pn- The validity of the quadratic approximation has 
been verihed with good numerical accuracy from isospin-symmetric nuclear matter [6 = 0) up to 
pure neutron matter (J = 1) by most of the existing nuclear many-body theories using various 
interactions [2]. Nonetheless, it has been shown consistently in numerous studies [3] that for some 
properties of neutron stars, such as the proton fraction at beta-equilibrium, the core-crust transition 
density and the critical density for the direct URCA process to occur, even a very small quartic 
isospin asymmetry energy A 4 (p) (multiplied with in the expansion of the energy per nucleon) 
can make a big difference. 

Given the fact that all the available numerical solutions of the nuclear many-body problem 
conhrm the validity of the quadratic approximation, the quartic isospin asymmetry A 4 (p) should 
be rather small. However, in the recent work by Cai and Li [1], which employs an empirically 
constrained isospin-dependent single-nucleon momentum-distribution and the equation of state of 
pure neutron matter near the unitary limit, a signiheant quartic isospin asymmetry energy of 
A/^{pq) = (7.2 ± 2.5) MeV has been found. This value amounts to about 16 times the free Fermi 
gas prediction, see eq.(2). On the other hand, the calculations by the Darmstadt group |5] based 
on chiral low-momentum interactions and many-body perturbation theory lead to a small value of 
A^{Pq) = (1.0 ±0.2) MeV. 

The purpose of the present paper it to give a prediction for the density-dependent quartic 
isospin asymmetry energy A 4 ^{kf) in the chiral approach to nuclear matter developed in refs.[6l [7]. 

^This work has been supported in part by DFG and NSFC (CRC 110). 
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In this approach the long- and medium-range NN-interactions arising from multi-pion exchange are 
treated explicitly and few parameters encoding the relevant short-distance dynamics are adjusted to 
bulk properties of nuclear matter. Single-particle potentials [7], quasi-particle interactions [S], the 
thermodynamic behavior of nuclear matter at hnite temperatures |H] and the density-dependence 
of the in-medium quark condensate [lO] follow then as predictions in that framework (see also the 
recent review article m)- 

The present paper is organized as follows. In section 2, analytical expressions are given for the 
contributions to the quartic isospin asymmetry energy A^lkf) as they arise from Ivr-exchange and 
chiral 27r-exchange. The three-nucleon interaction generated by 27r-exchange and excitation of a 
virtual A(1232)-isobar is considered as well. These interaction contributions lead at saturation den¬ 
sity po = 0.16 fm“^ (or kfo = 263 MeV) to the (small) value A 4 {kfo) = 1.5 MeV, which amounts to 
about three times the kinetic energy part. Moreover, in the course of the calculation one encounters 
components of the second-order Ivr-exchange whose representation of the fourth derivative with the 
respect to 5 at 5 = 0 is singular. In section3, the generic presence of a non-analytical term ^'^In |5| 
in the expansion of the energy per particle of isospin-asymmetric nuclear matter is demonstrated by 
calculating in closed form the second-order contribution from a s-wave contact-interaction. Clearly, 
after having established its existence, the non-analytical term In |(5| should be included in future 
hts of the equation of state of (zero-temperature) isospin-asymmetric nuclear matter. 


2 One-pion and two-pion exchange contributions 


In this section we collect the expressions for the quartic isospin asymmetry A 4 {kf) as they arise from 
one-pion and two-pion exchange diagrams following refs. PEI. Isospin-asymmetric (spin-saturated) 
nuclear matter is characterized by different proton and neutron Fermi momenta, kp^n = fc/(lT<^)^^^- 
Expanding the energy per particle at hxed nucleon density p = 2fcj/37r^ in the isospin asymmetry 
parameter 5 up to fourth order gives: 

E^{kp, kn) = E{kf) + S^A^ikf) + S^A^ikf) + 0{S^) , (1) 


with A 2 {kf) the (usual) quadratic isospin asymmetry energy. We view the density-dependent ex¬ 
pansion coefficients E{kf), A 2 {kf) and A 4 {kf) as functions of the Fermi momentum /cj, since they 
emerge in this form directly from the calculation. The hrst contribution to A 4 {kf) comes from the 
relativistically improved kinetic energy Tkin(p) = p^/2M — and it reads: 


A4{kfY^-^ = 


k) ( k) \ 
162MW^ ’ 


( 2 ) 


with M = 939 MeV the average nucleon mass. The corresponding value at nuclear matter saturation 
density po = 0.16 fm“^ (or fc/o = 263 MeV) is ^ 4 (^/ 0 )^'^™^ = 0.464 MeV. 

For the treatment of two-body interactions that depend on the momentum transfer |pi—P 2 I the 
following expansion formulas for integrals over two Fermi spheres are most helpful: 
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The z-dependent weighting functions at order and 6^^ have been obtained by applying several 
partial integrations. The contribution of the Ivr-exchange Fock diagram to the quartic isospin 
asymmetry energy reads: 
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with the dimensionless variable u = kf/m.,^. The second line in eq.(5) gives the relativistic 1/M^- 
correction. The occurring physical parameters are: nucleon axial-vector coupling constant qa = 1-3, 
(neutral) pion mass m-,, = 135 MeV, and pion decay constant /^r = 92.4 MeV. 

Next in the chiral expansion comes the iterated (second-order) Ivr-exchange. With two medium 
insertions |(1 -|- Tz)6{kp — |pj|) -|- |(1 — T^)6{kn — |pi|) one gets a Hartree contribution of the form: 
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and the corresponding Fock exchange-term reads: 
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Pauli-blocking effects at second order are included through diagrams with three (isospin-asymmetric) 
medium insertions |6]. We consider here only the factorizable Fock contribution for which the energy 
denominator gets canceled by factors from the momentum-dependent vriV-vertices (see eqs. (11,26) 
in ref. [6]). Its contribution to the quartic isospin asymmetry energy can be represented as a one- 
parameter integral, Ai^{k = 5f^Mm^(127r/^)“^ jQdxI{x, m), where the lengthy integrand I{x, u) 
involves the function ln[l + (m + a:)^] — ln[l + (m — a:)^] and its square. The corresponding value 
at saturation density is A 4 {kfQY^^^'> = —1.35 MeV, thus counterbalancing most of the Fock term 
^ 4 (^/ 0 )*'^^^ = 1.70MeV without Pauli-blocking written in eq.(7). For the non-factorizable pieces 
the representation of the fourth derivative with respect to 5 at 5 = 0 includes singularities of the 
form (m — x)~^,i> = 1,2. When subtracting these singular terms from the integrand only very 
small numerical values are obtained for the non-factorizable Hartree contribution. In the case 
of the quadratic isospin asymmetry energy A 2 {kfo) one hnds that the non-factorizable pieces (see 
eqs.(24,26) in ref. [6]) tend to cancel each other almost completely as: (—11.6 + 12.0) MeV. Therefore 
one can expect that the omission of the non-factorizable pieces does not change much the hnal result 
for the quartic isospin asymmetry energy A^lkf). However, the observation that the iterated Ivr- 
exchange has components with a singular representation of their fourth derivative with respect to S 
at 5 = 0, indicates that the expansion in eq.(l) becomes non-analytic beyond the quadratic order 
This feature is demonstrated in section 3 by calculating in closed form the second-order contribution 
from a s-wave contact-interaction. 

We continue with the contribution of the irreducible 2vr-exchange to the quartic isospin asym¬ 
metry energy. Using a twice-subtracted dispersion relation for the 2vr-exchange NN-potential in 
momentum-space and the master formulas in eqs. (3,4), one obtains: 
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Figure 1: Quartic isospin asymmetry A^lkf) as a function of the nucleon density p = 2 /cj/ 37 r^. 
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where ImVc,r and ImhFc; 7 ^ are the spectral functions of the isoscalar and isovector central and tensor 
NN-amplitudes, respectively. These imaginary parts are composed of the functions \J— 4m2 and 
arctan(y/i2 — 4m2/2A), with A = 293 MeV the delta-nucleon mass splitting. Note that due to 
the implemented subtractions the fcj-expansion of A^ik in eq.( 8 ) starts with the power fcj. A 
short-distance contribution proportional to kj is supplemented by the subtraction constants: 




Wkj / 2 B 5 
(3M)4 \Y' 


(9) 


with the parameters = 0 and Bn ^5 = —3.58 adjusted in ref. [7] to the empirical nuclear matter 
saturation point and quadratic isospin asymmetry energy = 34 MeV. 

Finally, we consider the long-range three-nucleon interaction generated by 27r-exchange and 
virtual excitation of a A(1232)-isobar [7]. The corresponding three-body Hartree contribution reads: 
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while the associated three-body Fock term can be represented as (127r/^)“"^A“^u), 
where the lengthy integrand J(x, u) involves the functions arctan(M -|- x) -|- arctan(M — x) and ln[l + 
{u -|- — ln[l + {u — xY]. Note that the three-body contact-term proportional to C, introduced 

additionally in ref. [7] does not contribute to the quartic isospin asymmetry energy Ai{kf). 

Summing up all the calculated contributions, one obtains the result for the density-dependent 
quartic isospin asymmetry energy A,^{kf) of nuclear matter as shown in Fig. 1 in the density region 
0 < p < 2po = 0.32 fm“^. The predicted value at saturation density po = 0.16 fm“^ is A 4 (fc/o) = 
1.49MeV. It amounts to 3.2 times the free Fermi-gas part = 0.464MeV. Note that 
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interaction contributions to A^i^kf) start (at least) with the power k^. The density-dependence of 
the full line in Fig. 1 is to a good approximation Actually, it should be noted that the present 
calculation of the quartic isospin asymmetry energy A 4 ^{kf) is performed in a framework where 
empirical constraints from bulk properties of nuclear matter are satished. This does not apply to 
the recent work in ref.[l], where a large value of A^lkfo) = (7.2 ± 2.5) MeV has been found. 


3 S-wave contact interaction to second order 


The analysis of the Pauli-blocking corrections to the second-order (iterated) l7r-exchange has indi¬ 
cated that non-analytical terms may occur in the ^-expansion of the energy per particle of isospin- 
asymmetric nuclear matter beyond the quadratic order. In the extreme case there could be a cubic 
term |5|^, which is after all even under the exchange of protons and neutrons: 6 —S. In order to 

clarify the situation, we consider a s-wave contact interaction: 


TT r 


Vet — [Os + Sot + {at — ttg) (Ti ■ a2 


( 11 ) 


and examine it in second-order many-body perturbation theory. For this simple interaction, the 
occurring integrals over three Fermi spheres with (at most) two different radii fcp, kn can be solved 
in closed analytical form. The pertinent function to express the result in the isospin-asymmetric 
conhguration of interest is: 


35 / dz {z — z‘^)\2xz -\- {x'^ — z^) In 1 

Jo I \x — z\) 

= -(15 + 33a;2 - Ax^) + \{A2x'^ - 15 - 35a;^) In + 2a;^ In —-r, (12) 

2 4 \x -1\ \x^ -1\ 


where the variable a: > 0 is set to a ratio of Fermi momenta, [(1 -f 5)/(l — 5)]^^/^ or 1. Note that the 
function dehned in eq.(12) has at a; = 1 the value 22 — 4 In 2. Combining the second-order Hartree 
and Fock diagrams generated by Vet according to their spin- and isospin-factors and performing the 
expansion in powers of 6, one obtains the following result for the energy per particle: 


F, (k k i(2nd) _ 


+ a?)(11 - 2ln2) + — [ 0^(3 - ln2) - 0^(2 + ln2) 


[a4l01n^ + 21n2-—) +a: 

8u 4 3 6 y 


301n^ + 21n2 + - 
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+ 0(^®) .(13) 


The crucial and novel feature which becomes evident from this expression is the presence of the non- 
analytical logarithmic term 5^1n(|5|/3). Interestingly, the corresponding coefficient is three times 
as large in the spin-triplet channel as in the spin-singlet channel. For comparison the first-order 
contribution of the s-wave contact interaction Vet reads, E^sikp, = /cj[ — as — at + 5‘^{at — 

as/3)]/27rM, without any higher powers of 6. Note that the sign-convention for the scattering 
lengths ag^t is chosen here such that positive values correspond to attraction. As a check we have 
rederived the same results at first and second order by using the alternative (and equivalent) form 
of the s-wave contact interaction V^, = 7r[3as + at + {as — at) ti ■ T 2 ]/M. 

In Fig. 2 the dependence of the second-order energy per particle Eas{kp^knY‘^°^"^ on the isospin 
asymmetry parameter 5 is shown for three different choices of the s-wave scattering lengths: ag = at, 
at = 0 and = 0. In each case the full line shows the exact result and the (nearby) dashed line 
gives the expansion in powers of 6 truncated at fourth order according to eq.(13). One observes 
that these expansions reproduce the full ^-dependence very well over the whole range — 1 < 5 < 1. 
Note also that the prefactor kjal^/bir'^M of dimension energy has been scaled out in Fig. 2. 
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Figure 2: Dependence of the 2nd order energy per particle Ei^s{kp, on the isospin asymmetry 

6. Three different choices for the scattering lengths, a* = at, at = 0 and = 0, are considered. 

If one performs for the second-order energy density the fourth derivative with respect to 5 at 5 = 0 
under the integral, then one encounters integrands with singularities of the form (1 — u = 1,2. 
The origin of these singularities, or in the proper treatment the non-analytical term <5^ ln(|(5|/3), lies 
in the energy denominator of second-order diagrams. For an inhnite (normal) many-fermion system 
the energy spectrum has a vanishing gap between bound states in the Fermi sea and excited states in 
the continuum. Such a gap-less energy spectrum causes a singularity, respectively a non-analyticity, 
if small asymmetries of the Fermi levels of two components are analyzed with too high resolution. 

In summary, we have demonstrated that the non-analytical term (5^1n(|5|/3) will be generically 
present in calculations of isospin-asymmetric nuclear matter that go beyond the mean-held Hartree- 
Fock level. Therefore, a term 5‘^ln|(5| should be included in future hts of the equation of state of 
(zero-temperature) isospin-asymmetric nuclear matter and its role should be further examined. 
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